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ABSTRACT 

Thermistor  measurements  of  small  scale  temporal  and  spatial  fluctu- 
ations of  temperature  were  made  near  the  surface  of  the  sea  from  the 
NUC  Oceanographic  Research  Tower.   Wave  height,  water  particle  velocity, 
sound  velocity,  and  salinity  were  measured  simultaneously.   Data  were 
subjected  to  statistical  and  spectral  analyses.   Temperature  fluctua- 
tions were  found  to  vary  with  wave  height.   A  linear  relation  was  found 
to  exist  between  temperature  fluctuations  and  wave  swell.   The  assumed 
relationship  was  verified  by  the  results  of  spectral  analysis.   High 
coherence  (.68)  between  wave  height  and  temperature  fluctuations 
occurred  near  the  frequency  of  maximum  wave  energy.   A  relatively  con- 
stant phase  difference  of  180  was  observed  in  the  frequency  band 
normally  associated  with  surface  swell  (.03  Hz  to  0.2  Hz).   Water 
particle  velocities  displayed  the  same  relatively  high  coherence  in  the 
frequency  band  of  high  wave  energy.   Temperature  fluctuations  were  180 
out  of  phase  with  the  u-component  of  velocity  and  90   out  of  phase  with 
the  w-component  of  velocity.   These  phase  differences  agreed  with  linear 
wave  theory.   Coherence  was  observed  to  be  lower  and  phase  difference 
more  random  with  increased  thermal  gradients.   This  was  attributed  to 
an  increase  in  turbulent  temperature  fluctuations. 
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I.   INTRODUCTION 

An  experiment  designed  to  measure  fluctuations  of  environmental 
parameters  in  a  small  volume  of  water  near  the  surface  of  the  sea  v?as 
conducted  on  8  and  9  June,  1972.   The  experiment  was  carried  out  on  the 
Naval  Undersea  Research  and  Development  Center's  Oceanographic  Research 
Tower  in  San  Diego,  California  by  members  of  the  Oceanography  Department, 
U.  S.  Naval  Postgraduate  School. 

The  main  objectives  of  the  portion  of  the  experiment  described 
herein  were  to: 

1.  Measure  the  spatial  and  temporal  temperature  field, 

2.  Define  a  relationship  that  describes  the  temperature  field 
in  the  near  surface  ocean, 

3.  Examine  the  relationships  between  the  temperature  field  and 
other  measured  parameters  via  spectral  and  statistical  analyses. 

Temperature  in  the  ocean  environment  may  be  regarded  as  a  passive 
scalar.   Passive  scalars  are  quantities,  including  contaminants,  that 
modify  and/or  can  be  transported  in  a  fluid  but  which  do  not  introduce 
bouyancy  effects.   A  scalar  which  introduced  bouyancy  effects  would  be 
considered  to  be  active.   Clearly  large  temperature  differences  can 
result  in  bouyant  motion  but  the  scales  of  temperature  fluctuations 
considered  here  render  such  induced  motion  negligible. 

Since  temperature  is  a  passive  scalar,  the  temperature  field  does 
not  create  motion  but  responds  to  motion.   The  temperature  field  is 
advected  and  convulsed  by  the  vector  velocity  field.   Hie  temperature 
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field  is  then  expected  to  reflect  changes  and  motion  occurring  within 
the  velocity  field.   As  such  it  is  a  good  indicator  of  the  types  of 
motion  occurring  in  the  near  surface  ocean. 
Temperature  change  can  be  caused  by: 

1.  the  vertical  motion  of  the  temperature  gradient  with  negligible 
changes  to  the  gradient, 

2.  the  movement  of  a  continuum  mass  of  water  with  distinctive 
temperature  structure, 

3.  wave  induced  motion, 

4.  turbulence. 

The  first  of  these  usually  result  from  internal  waves,  which  can 
occur  whenever  a  density  gradient  exists  in  the  fluid.   Internal  wave 
motion  can  be  irrotational ,  linear,  and  not  alter  the  mean  temperature 
gradient.   However  non-linear  internal  wave  motions  do  exist  and  can 
induce  large  changes  in  the  temperature  gradient.   A  large  body  of 
literature  on  internal  waves  exists  and  it  is  not  intended  to  discuss 
the  manifold  variety  of  these  waves  in  this  paper.   Only  motions  pertinent 
to  this  experiment  will  be  examined  in  order  to  interpret  results  of 
the  experiment. 

The  second  cause  of  temperature  fluctuations  is  due  to  water  with 
a  different  temperature  field  being  advected  past  the  point  of  interest, 
which  would  include  such  things  as  currents  and  large  eddies.   The 
contributions  from  this  source  are  considered  to  be  negligible  in  the 
area  of  this  experiment. 

Wave  induced  temperature  fluctuations  are  limited  to  depths  where 
wave  motion  is  in  evidence  which  normally  encompasses  depths  of  less 
than  one  wave  length.   Furthermore,  in  order  for  wave  action  to  cause 


temperature  fluctuations,  temperature  inhomogeneities  must  be  present  in 
the  water  which  surrounds  a  fixed  temperature  sensor. 

Oceanic  turbulence  is  a  broad  band  phenomenon  covering  scales  of. 
kilometers  to  fractions  of  a  centimeter.   It  is  thus  present  to  modify 
temperature  fluctuations  caused  by  each  of  the  above  three  processes. 
At  scales  sufficiently  removed  from  those  of  other  processes,  the 
temperature  field  itself  can  be  considered  to  be  a  turbulent  field  which, 
analogous  to  the  velocity  field,  has  an  "energy"  containing  region  and 
an  inertial  sub-range.   Due  to  the  difference  in  the  rates  of  viscous 
and  thermal  dissipation  in  the  ocean,  the  very  small  scale  temperature 
field  tends  to  persist  longer  than  the  velocity  field  and,  thus,  water 
which  was  at  one  time  turbulent  can  be  "marked"  by  a  turbulent  temper- 
ature field.   Although  such  small  scales  are  beyond  the  scope  of  the 
apparatus  used  for  this  experiment,  turbulence  at  larger  scales  will  be 
examined. 

In  order  to  meet  the  objectives  proposed  by  this  thesis,  it  is 
necessary  to  be  able  to  separate  the  various  causes  of  temperature 
fluctuations.   Fortunately  this  is  not  difficult  in  the  case  of  internal 
and  surface  waves.   These  tend  to  have  narrow  spectra  and  appreciable 
energy.   Furthermore  they  are  well  separated  from  each  other.   At  the 
experimental  site,  internal  waves  tend  to  be  present  at  periods  from 
4-10  minutes  whereas  a  typical  surface  wave  period  is  several  seconds. 
Resolution  of  advected  temperature  fluctuations  is  more  difficult  as 
many  scales  are  involved.   A  prior  model  of  the  temperature  field  must 
be  assumed  and  several  sensors  are  needed  to  adequately  define  horizontal 
and  vertical  scales. 
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The  time  scale,  length  scale,  measuring  device  used,  and  environ- 
mental conditions  all  effect  the  components  that  will  be  present  in  the 
measured  Eulerian  temperature  signal.   The  length  of  the  longest  data 
record  in  this  experiment  was  about  25  minutes,  the  spatial  separation 
of  sensors  ranged  from  23  centimeters  to  89  centimeters  and  sea  surface 
conditions  were  restricted  to  swell-type  waves.   The  above  conditions 
thus  restricted  the  main  causes  of  temperature  fluctuations  to: 

1.  Internal  waves 

2.  Surface   swell 

3.  Small    scale   turbulence 

Internal  waves  are  gravity  waves  that  can  exist  whenever  a  density 
gradient  exists  in  a  fluid.   In  the  ocean,  these  waves  propagate  along 
thermal  gradient  layers  where  the  vertical  change  in  density  is  very 
small  compared  to  the  average  density  in  the  water  column.   Under  these 
conditions  it  is  possible  to  have  very  long  internal  waves  with  large 
amplitudes  propagate  with  only  a  small  amount  of  energy.   Internal 
waves  are  most  often  manifested  in  oceanographic  data  as  Eulerian 
temperature  fluctuations.   Typical  internal  waves  for  the  area  of  this 
experiment  have  been ' described  in  the  literature  [Lafond , (7 ) ] .   He  found 
periods  ranging  from  four  to  eight  minutes,  amplitudes  from  60  to  240  cm., 
and  an  average  velocity  of  about  15  cm. /sec.   The  temperature  fluctu- 
ations associated  with  internal  waves  should  exhibit  a  linear  relation- 
ship with  the  instantaneous  wave  amplitude  of  the  internal  wave  at  the 
depth  under  consideration  and  the  temperature  gradient.   If  A  (t)  is  the 
instantaneous  amplitude  of  the  internal  wave  at  the  depth  under  consid- 
eration, Z  ,  then  a  reasonable  expression  for  temperature  fluctuations 
o 

due  to  internal  wave  propagation  is: 
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V  ■  a;  '  Ao(t)  <M> 


Surface  swell  is  known  to  be  closely  approximated  by  linear  theory. 
The  waves  for  this  study  were  analyzed  [Krapohl, (6) ]  and  found  to  obey 
linear  theory  to  a  high  degree  of  accuracy  and  thus  can  be  considered 
to  be  linear  and  irrotational .   If  Tl(t)  is  the  instantaneous  wave  height, 
then  the  temperature  fluctuations  due  to  surface  swell  should  be  of  a 
form  similar  to  those  due  to  internal  waves: 


Tj/t)   =  K(z,f)  '   g   *  n<t)  (1-2) 


where  K(z,f)  is  a  constant  that  relates  the  motion  at  the  surface  to 
the  motion  of  the  water  column  at  the  depth  under  consideration.   Surface 
swell  has  periods  ranging  from  6  to  20  seconds  and  amplitudes  on  the 
order  of  meters.   Furthermore,  from  equations  (1-1)  and  (1-2),  it  is 
clear  that  in  the  presence  of  a  negative  temperature  gradient  the  tem- 
perature fluctuations  will  be  180   out  of  phase  with  the  instantaneous 
amplitude  of  internal  and  surface  waves. 

The  separation  of  the  turbulent  temperature  signal  from  the  total 
signal  at  first  appears  quite  difficult.   Certainly  deterministic 
approaches  have  not  been  successful  in  turbulence  studies  and  statistical 
treatments  are  often  followed.   However  because  of  the  narrowness  of 
internal  and  surface  wave  spectra  the  dominating  cause  of  temperature 
fluctuations  is  usually  recognizable.   Furthermore  the  "energy"  contained 
in  turbulent  temperature  fluctuations  is  quite  small  which  means  that  at 
frequencies  where  internal  or  surface  waves  are  found,  turbulence  is  not 
a  major  feature. 
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In  order  to  determine  clearly  the  relationships  between  the  various 
parameters  measured  during  this  experiment,  it  is  desirable  to  remove 
the  effects  of  turbulence,  even  if  they  are  small,  from  the  signals  of 
the  other  parameters.   Spectral  analysis  partially  aids  vis  in  this 
direction  but  additional  analysis  schemes  need  to  be  investigated. 
Necessarily  such  schemes  will  be  band  limited  to  cover  the  scales  the 
experiment  is  able  to  resolve.   For  this  experiment,  the  turbulent 
frequency  range  would  be  from  about  0.5  Hz  on  up  to  the  Nyquist  frequency, 

Figure  (1)  displays  a  qualitative  prediction  of  the  thermal  energy 
density  distribution  that  could  reasonably  be  expected  under  the  con- 
straints of  this  experiment.   As  has  been  discussed,  where  surface 
swell  and  internal  waves  predominate,  there  are  also  "masked"  contribu- 
tions from  turbulent  processes. 
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Figure    1.      Predicted    energy   density   distribution   of 
temperature    fluctuations. 
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II.   THEORETICAL  CONSIDERATIONS 

Discussions  of  random  data  analysis  schemes  are  found  in  many 
references.   The  following  discussion  is  based,  in  large  part,  on 
chapters  3,5,6  and  9  of  Bendat  and  Piersol's  Random  Data,  Analysis  and 
Measurement  Procedures.   Throughout  this  discussion  temperature  fluctu- 
ations due  to  all  causes  are  assumed  to  be  stationary  and  independent 
from  one  another. 

Temperature  fluctuations  due  to  surface  swell  can  be  described  by 
a  constant  parameter  linear  system  whose  fundamental  properties  are 
invariant  with  time.   The  dynamic  characteristics  of  this  type  of  system 
can  best  be  described  through  the  use  of  a  weighting  function,  h(T) , 
which  is  defined  as  the  output  of  a  system  at  any  time  to  a  unit  impulse 
applied  at  a  time  T   before.   The  output,  T(t),  of  this  system  whose 
input  is  surface  wave  height,  T|(t),  is  given  by  the  convolution  integral: 

00 

T(t)   =  /*h(T)  •  (t  -  T)dT  (II-l) 

The  Fourier  transform  of  the  weighting  function  is  the  transfer 
function  defined  by: 

oo 

H(f)      =      f  h(T)e~i2rrftdT    ;  (II-2) 

''o 

where    f   is    frequency.      From   equation    (II-l)    it   can   easily  be   shown   that: 

oo 

T(t   +  T)    =      /h(T)11(t   +  T    -    OdC    '»  <II_3> 

o 

where  £  is  a  dummy  variable. 
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The    product   of   T(t)    and    T(t  +  T)    is    given  by: 

CO 

T(t)    •    T(t  +  T)    =ffh(C)T\(t   -  C)    '    h(e)T)(f  +  T    -   e)ded£  ;         (II-4) 

where  £  and  e  are  dummy  variables. 

The  time  average  of  the  above  equation  is  defined  as  the  temporal 
auto-covariance  function,  R™(7) .   In  terms  of  weighting  functions  and 
the  input  auto-covariance  function,  R_(T)  is  given  by: 

CO 

VT)  B//h<e>h<OYT  +  e  -Oded<  (II-5> 

The  cross-covariance  function  between  the  input  and  output  can 
similarly  be  derived  using  as  a  starting  point  the  relation: 

CO 

T](t)  ■  T(t  +  T)  -  fh(e)T\(t)    .  T](t  +  T  -  e)de  (II-6) 

o 

which  yields  on  time  averaging  the  cross-covariance  relation: 

CO 

V(T)   =  /h(e)YT  -  c>de  (II"7) 

The  Fourier  transform  of  the  auto-covariance  function  is  related  to 
the  spectral  energy-density  function  by  the  relation: 

CO 

ST(f)      =      bj      R^T)    cos(2rrfT)dT  (II-8) 
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The  Fourier  transform  of  the  cross-covariance  function  is  similarly 
related  to  the  cross-spectral  energy  density  function  by  the  relation: 


00 


v(f)  =  /v(r)e"i2nft  dT  (II'9) 

From  equations  (II-2) , (II-8) ,  and  (II-9)  two  important  equations 
in  the  frequency  domain  can  be  derived: 


sT(f)  - 


H(f) 


2  •  S^Cf)  (11-10) 


and 


S^f)  =  H(f)  •  S  (f)  (11-11) 


In  the  frequency  band  of  swell  waves  (about  .05  Hz  to  .3  Hz)  tem- 
perature fluctuations  are  composed  primarily  of  fluctuations  due  to 

surface  waves,  T  ,  and  fluctuations  due  to  turbulence,  T  .   Thus  the 
w  t 

temporal  temperature  field,  T(t) ,  can  be  linearly  expressed  as: 


T(t)   =  Tw(t)  +  Tfc(t)  (11-12) 


It  is  assumed  that  the  wave  induced  and  turbulent  temperature 
fluctuations  are  statistically  independent.   The  spectral  energy  density 
function  of  the  temperature  fluctuations  is  then  given  by: 


S  (f)   =   S   (f)   +   S   (f)  (11-13) 

W  t 
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The  coherence  function,  y  _(f),  between  wave  height  and  temperature 
is  defined  by: 


Vf)  = 


Vf> 


S  (f).sT(f) 


(11-14) 


Assuming  that  T   and  T   are  statistically  independent  implies  that 

W         T- 


s^T(f) 


\ 


(f) 


(H-15) 


W 


If  it  is  further  assumed  that  the  temperature  fluctuations  due  to 
waves  can  be  expressed  as  a  constant  parameter  linear  system  with  wave 
height  as  the  input,  then  clearly  the  coherence  between  wave  height  and 
temperature  should  be  equal  to  one: 

2 

=   1  (11-16) 


S    (f) 

1  W 


yT]T  (f)       ST  (f)  •  S  (f) 
W  Tw       Tl 


Substituting  equations  (11-13),  (11-15)  and  (11-16)  into  equation 
(11-14)  and  simplifying  yields: 


1 


sT  (f) 

w 


(11-17) 


The  assumption  that  the  turbulence-  and  wave-  generated  temperature 
fluctuations  are  statistically  independent  and  uncorrelated  allows  us  to 
arrive  at  equation  (11-17) .   Allowing  correlation  between  the  two 
variables  would  add  to  the  complexity  but  would  still  yield  the  same 
basic  result  and  thus  the  simplification  is  preferred  for  discussion 
purposes. 
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From  equation  (11-17)  it  can  be  seen  that  turbulence-generated 
temperature  fluctuations,  which  can  be  considered  as  noise  in  the  pro- 
cess, cause  the  coherence  between  wave  height  and  temperature  fluctua- 
tions to  be  less  than  one.   Furthermore  it  is  also  obvious  that  if  the 
spectrum  of  wave  generated  temperature  fluctuations  is  held  nearly 
constant  then  an  increase  in  turbulence  will  cause  the  coherence 
between  wave  height  and  temperature  fluctuations  to  become  smaller. 

If  the  assumed  relationship  between  surface  waves  and  temperature 
fluctuations  is  correct,  i.e.; 


Tw(t)   =  K(z,f)  ■  ~   ■  T](t)  (11-18) 


then  we  would  expect  that  in  the  absence  of  turbulence  the  coherence 
function  would  reach  a  value  of  one  at  the  frequency  of  the  predominant 
swell  and  furthermore  that  in  the  presence  of  a  negative  vertical  tem- 
perature gradient  we  should  find  a  phase  angle  between  waves  and  tem- 
perature fluctuations  of  180  . 

Further,  it  should  be  noted  that  if  turbulent  temperature  fluctuations 
are  present  then  the  coherence  will  be  less  than  one  as  shown  by  equation 
(11-17).   Also  turbulence  would  be  expected  to  modify  the  phase  angle  to 
some  extent  in  a  more  or  less  random  fashion. 
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III.   INSTRUMENTATION 

The  temperature  sensing  device  chosen  for  this  experiment  was  the 
type  K496  double  beaded  thermistor  manufactured  by  Fenwall  Electronics 
Inc.,  Framingham,  Massachusetts.   This  thermistor  has  an  exceptionally 
thin  glass  coating  which  provides  a  water-tight  covering  and  allows  for 
a  low  thermal  response  time  of  less  than  150  milliseconds.   The  resistance 
of  the  thermistor  is  about  18  K-ohms  at  20   C  and  varies  about  870  ohms 
per  1   C  in  the  temperature  range  from  16   C  to  22   C.   These  high 
resistance  values  give  the  thermistor  a  high  degree  of  resolution.   The 
thermistor  also  has  very  low  self-heating  and  thus  in  seawater  can  pro- 
vide temperature  measurements  accurate  to  .01   C.   Further  details  on 
the  type  K496  can  be  found  in  reference  (10). 

The  thermistor  was  supplied  mounted  with  epoxy  potting  compound  in 
a  3mm  diameter  hollow  steel  tube.   The  leads  of  the  thermistor,  which 
extended  from  the  tube,  were  soldered  to  plastic-coated,  shielded 
instrumentation  cable.   This  connection  and  the  end  of  .the  steel  tube 
were  placed  in  a  mold  which  was  then  filled  with  epoxy  potting  compound. 
On  hardening,  the  epoxy  formed  a  6  cm  in  diameter  base  for  the  steel  tube 
which  protruded  about  9  cm  above  the  base.   The  overall  sensor  length 
was  about  14  cm.   A  total  of  six  sensors  were  prepared  and  calibrated 
in  a  temperature  controlled  bath  using  a  previously  calibrated  digital 
quartz  thermometer  over  the  range  from  10   C  to  22   C.   A  typical  tem- 
perature versus  resistance  calibration  curve  is  shown  in  Figure  (2). 
For  any  small  temperature  range  (±  2      C)  the  curve  is  linear  to  a  high 
degree  of  accuracy. 
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Each  sensor  was  used  as  one  leg  of  a  Wheatstone  bridge  (Figure  3). 
The  bridge  consisted  of  two  24.3  K-ohm  resistors,  the  thermistor,  and 
a  variable  resistor  which  was  used  to  balance  the  bridge  to  zero  voltage 
output  at  any  selected  temperature.   The  output  of  the  bridge  was  thus 
a  direct  measurement  of  temperature  fluctuations  about  a  selected  mean. 
Typical  values  of  bridge  output  for  a  selected  mean  temperature  of  18   C 
were  305  mv  per  degree  C.   The  advantages  of  this  type  of  sensor  are 
high  sensitivity,  high  accuracy,  interchangeability ,  fast  thermal  response 
time  and  low  cost.   The  only  major  disadvantage  for  oceanographic  field 
work  is  that  they  are  very  fragile. 

In  addition,  wave  height,  water  particle  velocity,  sound  velocity, 
and  salinity  were  measured.   Wave  height  was  measured  using  the  Baylor 
wave  guage  which  has  an  accuracy  of  1  percent.   Water  particle  velocities 
were  measured  on  the  Engineering  Physics  Company's  Water  Current  Meter 
Model  EM CM  3B.   The  time  constant  of  this  device  is  0.2  seconds  and 
velocity  is  measured  to  1  percent  of  full  scale  deflection  which  for 
this  experiment  was  1  cm/sec.   These  sensors  are  described  in  detail 
by  Krapchl  [6],   Sound  velocity  was  measured  on  the  Ramsey  Corporation 
Mark  I  SVTD  [Gossner  (5)].   Salinity  was  measured  on  the  Bisset-Berman 
STD  [Frigge  (4)]. 
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IV.  FIELD  EXPERIMENT 

The  site  chosen  for  conducting  this  experiment  was  the  Naval 
Undersea  Research  and  Development  Center's  Oceanographic  Research  Tower. 
The  tower  is  located  in  18  meters  of  water  approximately  one  mile 
offshore  of  Mission  Beach,  San  Diego,  California.   The  tower  is  a  steel 
and  concrete  structure  mounted  on  four  tubular  steel  legs  which  are 
rigidly  attached  to  the  bottom.   An  enclosed  upper  deck  provides 
adequate  space  for  the  installation  of  a  wide  variety  of  electronic 
equipment  and  a  stable  power  supply  is  provided  via  underwater  cable 
from  shore.   Rail-mounted  carriages  are  mounted  on  the  north,  south, 
and  west  sides  of  the  tower  and  can  be  moved  up  and  down  the  sides  of 
the  tower  with  electric  winches.   Up  to  1000  pounds  of  sensors  can  be 
mounted  on  these  carriages  and  measurements  can  be  made  at  any  desired 
depth  from  the  surface  to  one  meter  above  the  sea  floor. 

An  instrument  mounting  frame  was  constructed  to  fit  the  carriages. 
Sensors  were  mounted  on  the  frame  with  consideration  given  to  minimizing 
interference  between  sensors  while  ensuring  good  coverage  of  spatial 
scales  expected.   Figure  (4)  shows  schematically  the  positioning  of  the 
various  sensors.   The  frame  and  sensors  were  mounted  on  the  west  side  of 
the  tower  and  were  thus  open  to  the  full  effects  of  the  swell  which  was 
propagating  from  west  to  east. 

Instrumentation  wire  was  led  from  each  sensor  to  associated  elec- 
tronics in  the  enclosed  upper  level.   The  output  of  the  Wheatstone 
bridge  associated  with  each  temperature  sensor  was  connected  to  a  Preston 
Model  8300  XWB  operational  amplifier  where  a  gain  factor  of  ten  was  applied 
prior  to   recording  the  signal  on  the  Sangamo  Model  3562,  14-channel 

magnetic  tape  recorder. 
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Figure  A.       Schematic   of   sensor   placement 
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On  8  and  9  June,  1972,  measurements  were  made  at  ten  preselected 
depths.   Outputs  from  each  sensor  were  recorded  on  eleven  of  the  four- 
teen available  tape  recorder  channels.   Major  difficulties  were  en- 
countered with  the  temperature  sensors.   Failure  of  three  of  the  four 
available  sensors  took  place  early  in  the  experiment.   Sensor  T4 
(Figure  4)  failed  at  the  beginning  of  the  first  run  and  sensors  T2  and 
T3  failed  early  in  the  second  run.   In  addition,  the  record  from  sensor 
T2  was  interupted  for  approximately  one  minute  in  the  middle  of  the 
first  run  which  effectively  removed  it  from  any  statistically  accurate 
analysis  scheme.   Thus  no  spatial  information  could  be  resolved  for 
distances  less  than  89  centimeters,  the  distance  between  sensors  Tl  and 
T3. 

The  thermistor  beads  of  sensors  T2,  T3  and  T4  were  found  to  have 
sustained  major  damage  and  the  most  probable  cause  of  this  damage  was 
contact  with  the  heavy  cable  of  the  Baylor  wave  guage.   In  normal 
operation  this  cable  is  placed  in  tension  but  in  order  to  keep  the  wave 
guage  centered  over  the  sensor  mounting  rack  it  was  necessary  to  attach 
one  end  of  the  cable  to  the  rack  and  thus  the  heavy  cable  had  to  be 
raised  and  lowered  by  hand  through  the  upper  few  meters  of  water  prior 
to  tension  being  applied.   It  is  significant  that  the  only  sensor  to 
remain  in  operation  for  the  duration  of  the  experiment  was  located 
farthest  from  the  tower  and  in  a  position  where  it  was  unlikely  to 
make  contact  with  the  wave  guage. 
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V.   DATA  CONVERSION  AND  ANALYSTS 

The  analog  data  record  was  played  back  through  a  Ci  5000  computer 
onto  an  eight  channel  strip  chart  recording  for  editing  purposes.   Strip 
chart  recordings  of  each  run  were  first  checked  for  errors  that  would 
make  them  unsuitable  for  digital  analyses.   These  included  such  things 
as  short  records,  clipped  signals,  and  obvious  electronic  malfunctions. 
The  runs  that  were  not  rejected  by  this  procedure  were  then  subjectively 
analyzed  on  the  basis  of  environmental  parameters,  presence  or  absence 
of  internal  wave  induced  temperature  fluctuations,  time  of  data  record- 
ing and  sensor  depth.   Four  runs  were  subsequently  chosen  for  analysis 
and  are  summarized  in  Table  I.   The  length  given  for  each  record  is  the 
length  of  the  continuous  record  used  for  digital  analysis  after  each 
signal  was  edited  to  ensure  that  all  systems  had  stablized.   Environ- 
mental conditions  that  prevailed  during  each  run  are  summarized  in 
Table  II.   In  addition  to  the  parameters  listed,  vertical  temperature 
structure  for  each  run  was  measured  with  a  mechanical  bathythermograph 
at  the  start  of  each  run  (Figure  5).   No  values  are  placed  on  the  tem- 
perature axes  of  this  figure  as  the  mechanical  BT  does  not  yield 
accurate  absolute  temperatures.   However  it  does  yield  accurate  tem- 
perature gradient  with  depth  information.   The  data  show  that  from 
Run  1  (1405  Local  Time)  to  Run  6  (1922  Local  Time)  the  vertical  tem- 
perature gradient  was  reduced  to  an  isothermal  condition  near  the  sur- 
face.  Winds  of  from  7-10  knots  from  the  southwest  were  recorded  from 
several  hours' before  Run  1  through  Run  5.   This  wind  speed  is  sufficient 
to  cause  some  mechanical  mixing  of  the  surface  waters.   Also  the  air 
temperature  decreased  from  67.5   F  to  64   F  from  Run  1  to  Run  6.   This 
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TABLE  I   Summary  of  Runs  Selected  for  Analysis 


Run  # 


Local 

Time/Date 


Length 
(Min.) 


Depth  of 
Sensors  (Ft.) 


1405/8  June 
1830/8  June 
1922/8  June 
0837/9  June 


19.25 
16.67 
17.81 
27.77 


18.7 
21.8 
14.3 
14.3 


TABLE   II      Summary   of    Environmental   Conditions    for    Selected    Runs 


Cloud 

Wind 

w 

ind  Spe 

ed 

Air  Temp. 

Sea  Sfc. 

Run  # 

Cover  (7o) 

Dir. 

(KTS) 

(F) 

Temp.  (F) 

1 

100 

WSW 

8 

67.5 

66.5 

5 

100 

SW 

7 

64.2 

66.0 

6 

100   • 

SW 

4 

64.0 

65.7 

8 

100 

SW 

4 

63.1 

65.7 
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Figure  5.   Bathythermographs  for  selected  runs. 
( indicates  sensor  depth) 


29 


decrease  was  accompanied  by  a  decrease  in  sea  surface  temperature  from 
66.5   F  to  65.7   F.   This  surface  water  cooling  would  be  expected  to  be 
accompanied  by  convective  overturn  and  therefore  mixing  of  the  surface 
waters.   Either  of  these  processes  would  give  rise  to  the  observed 
vertical  thermal  structure  changes  and  most  likely  both  were  at  work. 
Both  processes  are  non-linear,  turbulent  processes.   From  the  beginning 
of  Run  6  till  the  end  of  Run  8  winds  were  calm.   The  increased  near- 
surface  vertical  thermal  gradient  observed  in  Run  8  is  probably  due  to 
warming  of  the  surface  waters  by  the  early  morning  sun,  although  the 
sky  remained  overcast  throughout  the  experiment. 

Each  of  the  runs  selected  for  analysis  were  digitized  by  inputing 
the  analog  record  into  a  hybrid  system  consisting  of  a  Ci  5000  analog 
computer  and  a  XDS  9300  computer.   A  standard  analog  to  digital  con- 
version program  was  used  to  generate  a  7-track  digital  magnetic  tape. 
Since  a  great  deal  of  core  space  was  required  for  digital  analysis  of 
the  data  the  7-track  tape  was  converted  to  a  9-track  tape  which  was 
compatible  with  the  IBM  360/67  computer.   Conversion  was  executed  on 
the  IBM  360/67  using  a  standard  conversion  program.   Analysis  of  each 
selected  run  was  then  performed. 

For  each  pair  of  variables  (e.g.,  temperature  and  sound  velocity 
fluctuations) ,  means  and  variances  were  calculated  and  the  means  were 
subtracted.   Any  linear  trends  were  removed  and  auto-covariance  and 
cross-covariance  functions  were  calculated.   A  Parzen  window  was  applied 
to  the  covariance  functions  to  account  for  the  finite  length  of  records. 
The  spectral  energy  density  spectra  and  the  cross-spectral  energy  density 
spectrum  were  calculated  by  Fourier  transforming  the  auto-  and  cross- 
covariance  functions.   The  phase  and  coherence-functions  were  calculated 
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using  the  cross-spectral  quantities.   Graphs  of  the  auto-correlation, 
cross-correlation,  spectral  energy  density,  coherence,  and  phase 
functions  were  plotted  using  the  Cal-comp  plotter. 
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VI.   DISCUSSION  OF  RESULTS 

A.   SURFACE  WAVE  INDUCED  TEMPERATURE  FLUCTUATIONS 

Figure  (6)  is  a  plot  of  energy  density  spectra,  coherence  and  phase 
for  temperature  and  wave  height  for  Run  6.   Maximum  energy  of  the 
temperature  spectrum  is  centered  at  the  frequency  of  maximum  surface 
wave  energy  (.064  Hz)  which  corresponds  to  a  wave  period  of  approximately 
15.6  seconds.   Maximum  coherence  (.68)  is  found  at  this  frequency  and 
relative  maxima  in  coherence  generally  occur  at  frequencies  corresponding 
to  those  of  relative  maxima  in  the  wave  spectrum.   The  phase  relation- 
ship is  relatively  constant  in  the  range  from  about  .03  Hz  to  .2  Hz  at 
a  value  of  180  .   The  waves  in  this  frequency  band  were  swell  type  waves. 
The  high  coherence  values  at  frequencies  of  relative  maxima  in  the  wave 
spectrum  and  the  180   phase  difference  in  the  swell  frequency  band 
apparently  verify  the  assumed  relationship  between  surface  swell  and 
temperature  fluctuations.   The  coherence  and  phase  relationships  between 
temperature  fluctuations  and  water  particle  velocity  components  further 
verify  the  temperature/wave  relationship.   Linear  theory  predicts  that 
the  horizontal  components  of  wave-induced  water  particle  velocity  should 
be  in  phase  with  wave  height  fluctuations  and  that  the  vertical  velocity 
component  should  be  90   out  of  phase.   As  previously  stated,  Kraphol  [6] 
showed  that  the  waves  associated  with  this  experiment  were  highly  linear. 
The  phase  relations  outlined  above  should  hold  between  temperature 
fluctuations  and  the  velocity  components.   Figures  7  and  8  are  plots  of 
spectra,  coherence,  and  phase  for  T  vs  u  and  T  vs  w  respectively.   T  vs  u 
displays  relatively  high  coherence  (.61)  at  the  frequency  of  maximum 
wave  energy  with  a  phase  angle  of  180  .   T  vs  w  displays  high  coherence 

at  the  same  frequency  but  with  a  phase  angle  of  90°  from  .03  Hz  to  .22  Hz, 
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Figure   6. 

Spectra,    coherence,    and   phase 
for   temperature  and  wave  height, 
Run   6 
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Figure   7. 

Spectra,    coherence,    and   phase 
for    temperature   and   u-component, 
Run  6 
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Figure  8. 

Spectra,  coherence,  and  phase 
for  temperature  and  v-conponent , 
Run  6 
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An  anomaly  in  this  pattern  exists  at  about  .1  Hz.   This  frequency 
corresponds  to  a  relative  maximum  in  the  energy  density  spectrum  of 
temperature  fluctuations  and  a  relative  minimum  in  the  spectra  of  wave 
height,  u-  and  w-compononts  of  velocity.   A  coherence  maximum  (.2) 
occurs  at  this  frequency  for  both  T  vs  w  and  T  vs  wave  height.   The 
phase  between  temperature  and  wave  height  is  -120  which  does  not  agree 
with  the  predicted  phase  angle.   The  phase  angle  between  temperature  and 
w-component  of  velocity  is  150   and  thus  the  90   phase  angle  between 
wave  height  and  w-component  of  velocity  is  as  predicted  by  linear  theory. 

Figures  9  through  17  are  similar  plots  for  temperature,  wave  height, 
and  water  particle  velocities  for  Runs  5,  1  and  8.   All  runs  are  in 
general  agreement  with  the  results  of  Run  6.   Run  5  has  almost  identical 
results,  displaying  slightly  lower  coherence  values  at  the  frequency 
of  maximum  energy  and  slightly  less  consistency  in  the  phase  relation- 
ships.  These  results  still  show  very  good  adherence  to  the  assumed 
temperature-surface  wave  relationship.   Runs  1  and  8  show  significantly 
lower  coherence  values  and  more  random  phase  relationships.   These 
results,  though  less  striking,  are  still  in  general  agreement  with  those 
of  Runs  5  and  6. 

Figure  5  provides  a  possible  explanation  for  the  differences  observed 
between  Runs  5  and  6  and  Runs  1  and  8.   It  was  shown  in  equation  (11-14) 
that  in  the  frequency  band  of  surface  swell,  the  coherence  between  tem- 
perature and  wave  height  will  be  strictly  less  than  one  if  temperature 
fluctuations  are  assumed  to  be  composed  of  wave  induced  and  turbulent 
fluctuations.   Furthermore  it  was  shown  that  as  turbulent  fluctuations 
increased,  the  coherence  would  decrease.   Since  the  swell  waves  present 
remained  relatively  constant  for  all  runs,  we  would  expect  that  at 


36 


in 
z: 

LU 

cn 
> 

IS) 

cr 
z: 

LU 


cr. 
cr 

h- 

LJ 

LU 

CL 
CO 

L^ 

_J 


Figure   9. 

Spectra,  coherence,  and  phase 
for  temperature  and  wave  height, 
Run  5 
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Figure   10. 

Spectra,    coherence,    and   phase 
for   temperature  and  u-component, 
Run  5 
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Spectra,  coherence,  and  phase 
for  temperature  and  w-component, 
Run  5 
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Spectra,    coherence,    and   phase 
for    temperature  and  wave  height, 
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Spectra,    coherence,    and   phase 
for    temperature  and   u-component, 
Run   1 
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Spectra,    coherence,    and   phase 
for   temperature  and  v-component. 
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Spectra,  coherence,  a:id  phase 
for  temperature  and  vzave  height, 
Run  8 
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Spectra,    coherence,    and   phase 
for    temperature  and   u-component, 
Run   8 
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Spectra,    coherence,    and   phase 
for   temperature   and   v-component, 
Run  8 


Kid    ; 


zz 

cr   k 

LiJ     * 
X 

o 
u 


IK10 


FREQUENCY    (HEJ 


45 


constant  depth,  the  turbulent  velocity  field  strength  would  remain 
relatively  constant  in  the  frequency  band  under  consideration.   Since 
Run  6  and  Run  8  were  made  at  the  same  depth,  the  differences  between 
them  cannot  be  attributed  to  changes  in  the  turbulent  velocity  field. 
In  the  absence  cf  a  well  defined  temperature  gradient  (Runs  5  and  6)  we 
would  expect  turbulent  temperature  fluctuations  to  be  quite  small;  and 
with  a  larger  gradient  (Runs  1  and  8),  we  would  expect  higher  levels  of 
turbulent  temperature  fluctuations  regardless  of  the  turbulent  velocity 
field.   For  example,  we  would  observe  no  turbulent  temperature 
fluctuations  in  a  turbulently  well  mixed,  constant  temperature  bath. 

B.   INTERNAL  WAVE  INDUCED  TEMPERATURE  FLUCTUATIONS 

Inspection  of  strip  chart  recordings  showed  that  Run  1  was  the  only 
run  selected  for  analysis  that  contained  any  significant  long  period 
components.   The  temperature  signal  from  thermistor  Tl  was  summed  and 
averaged  for  each  12  second  interval  for  all  records  analyzed.   The 
results  of  the  averaged  temperature  signal  for  Run  1  are  plotted  in 
Figure  18.   Only  three  long  term  fluctuations  were  observed.   The 
fluctuations  displayed  a  6  to  9  minute  period  and  were  quasi-cycloidal 
in  shape.   The  average  temperature  change  from  maximum  temperature  to 
minimum  temperature  was  approximately  .48   C.   If  the  temperature 
fluctuations  due  to  internal  waves  do  in  fact  reflect  the  vertical 
movement  of  the  vertical  temperature  gradient,  [Equation  I-l],  then 
with  a  measured  thermal  gradient  of  .008   C/cm,  this  temperature 
fluctuation  corresponds  to  an  internal  wave  height  of  60  cm.   It  should 
be  noted  that  negative  temperature  fluctuations  indicate  positive  wave 
heights  and  vice-versa  if  the  assumed  internal  wave-temperature 
relationship  is  correct.   The  phase  difference  with  a  negative  gradient 
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Figure  18.   Temperature  averaged  each  12  seconds  for  Run  1 


47 


should  be  180  .   Figure  19  is  a  plot  of  the  v-velocity  component  averaged 
each  12  seconds  in  the  same  way  as  the  average  temperature  plot.   On  this, 
the  average  temperature  fluctuation  has  been  superimposed  and  the  phase 
reversed  180  .   The  same  internal  wave  can  be  seen  to  manifest  itself 
in  the  water  particle  velocity  record.   Referring  to  Figures  6  through 
17,  it  can  be  seen  that  for  Run  1  the  coherence  values  at  low  frequencies 
are  very  high  for  T  vs  u  (.72)  and  T  vs  v  (.62)  and  the  phase  angle  is 
180  .   For  Runs  5,  6  and  8,  which  had  little  internal  wave  activity, 
coherence  values  are  quite  low  at  low  frequencies.   Figures  20  and  21 
are  plots  of  spectra,  coherence  and  phase  for  sensor  T3  and  the  u-  and 
v-components  of  velocity  for  Run  1.   Values  of  coherence  and  phase  at 
low  frequencies  are  observed  to  be  similar  to  those  obtained  from  sensor 
Tl. 

C.   SPATIAL  RELATIONSHIPS 

Due  to  breakage  of  thermistors  and  technical  problems  with  electronics, 
only  Run  1  has  a  record  where  more  than  one  temperature  sensor  functioned. 
Temperature  sensors  Tl  and  T3  (Figure  4)  functioned  throughout  Run  1. 
The  distance  between  these  sensors  was  89  cm.   Figure  22  is  a  plot  of  the 
auto-correlation  functions  for  sensors  Tl  and  T3 .   They  are  quite  similar 
out  to  about  a  lag  of  22  seconds,  after  which,  there  is  a  difference 
between  them  of  2%  and  more,  but  they  are  in  general  agreement.   Figure  23 
is  a  plot  of  the  cross-correlation  between  Tl  and  T3 .   It  has  a  maximum 
value  of  .95  at  a  lag  time  of  14  seconds.   Thus  Tl  and  T3  are  most  highly 
correlated  when  lagged  14  seconds.   Figure  24  is  a  plot  of  the  spectra, 
coherence  and  phase  of  Tl  and  T3 .   The  highest  coherence  value  (.93)  is 
found  in  the  low  frequencies.   The  14  second  lag  for  a  maximum  cross- 
correlation  value  is  thus  probably  associated  with  internal  waves.   If 
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Figure  19.   V-component  of  velocity  averaged  each  12  seconds  with  averaged 
temperature  signal  superimposed  and  shifted  180  . 
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Figure   20. 

Spectra,    coherence,    and   phase 

for   temperature    (T3)    and   v-component, 

Run   1 
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Figure  21. 

Spectra,  coherence,  and  phase 

for  temperature  (T3)  and  u-component , 

Run   1 
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Fieure  22.   Auto-correlation  functions  for  thermistors  Tl  and  T3,  Run  1 
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Figure  24. 

Spectra,  coherence,  and  phase 
for  thermistors  Tl  and  T3, 
Run  1 
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we  assume  that  14  seconds  is  the  time  required  for  an  internal  wave  to 
travel  from  Tl  to  T3,  than  the  average  internal  wave  velocity  is  6.35 
cm/sec.  which  agrees  fairly  well  with  other  data  [Lafond  (7)].   A  second 
maximum  value  (.63)  occurs  in  the  coherence  between  Tl  and  T3  at  the  pre- 
dominant wave  frequency  of  .065  Hz.   The  energy  density  spectra  of  Tl 
and  T3  are  close  to  identical  to  a  frequency  of  about  .09  Hz.   The  fact 
that  the  coherence  is  less  than  one  is  attributed  to  turbulence. 

D.   RELATIONSHIP  BETWEEN  TEMPERATURE  FLUCTUATIONS  AND  SOUND  VELOCITY 

Sensor  T3  was  located  23  cm  from  the  sound  velocimeter.   Figure  25 
is  a  plot  of  spectra,  coherence  and  phase  for  T3  and  sound  velocity  for 
Run  1.   Fairly  high  values  of  coherence  are  apparent  to  a  frequency  of 
.12  Hz.   At  frequencies  below  about  .01  Hz,  coherence  is  close  to  1.0 
and  T3  and  sound  velocity  remain  in  phase  to  about  .02  Hz.   Near  the 
predominant  wave  frequency,  coherence  reaches  a  value  of  .72  and  phase 
is  relatively  constant  at  36   -  40  .   Little  coherence  is  exhibited 
above  .1  Hz.   Figure  26  is  a  similar  plot  for  Tl  and  sound  velocity. 
Tl  was  located  66  cm  farther  away  from  the  sound  velocimeter  than  T3 . 
In  general  the  coherence  between  Tl  and  sound  velocity  is  similar  to  that 
between  T3  and  sound  velocity  but  at  significantly  lower  levels.   Below 
.01  Hz  coherence  is  reduced  to  .92  and  at  the  frequency  of  predominant 
waves  coherence  is  .59.   In  addition,  the  band  of  frequencies  in  which 
temperature  and  sound  velocity  remain  coherent  is  much  narrower  for  Tl 
than  for  T3 .   The  differences  between  the  two  temperature  vs  sound 
velocity  coherence  plots  may  be  attributed  to  the  existence  of  turbulent 
fluctuations  between  the  two  sensors  or  the  fact  that  the  time  constants 
of  the  sound  velocimeter  and  the  thermistor  were  different. 
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Figure   25. 

Spectra,  coherence,  and  phase 
for  T3  and  sound  velocity, 
Run  1 
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Figure   26. 

Spectra,  coherence,  and  phase 
for  Tl  and  sound  velocity, 
Run  1 
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VII.  CONCLUSIONS 

1.   The  relationship  between  temperature  fluctuations  and  surface 
wave  height  (at  least  for  swell-type  waves)  is  verified  to  be  of  the 
form : 


T(t)   =  K(f,z).  •  g   .   T](t) 


2.  All  evidence  points  toward  a  similar  relationship  between 
internal  waves  and  temperature  fluctuations.   Only  one  twenty  minute 
record  displayed  sufficient  internal  wave  activity  to  verify  this. 

3.  Turbulent  temperature  fluctuations  greatly  degrade  the  coherence 
between  two  sensors  located  only  89  cm  apart. 

4.  Further  experimentation  is  needed  to  establish  the  relationship 
between  non-swell  surface  waves  and  temperature  fluctuations  and  to 
further  verify  the  internal  wave  relation. 

5.  To  successfully  measure  small  scale  spatial  correlations  of 
the  temperature  field,  sensors  must  be  placed  at  several  points  with 
separations  less  than  89  cm.   Furthermore  in  order  to  do  this  with  fast 
response  thermistors  a  great  deal  of  redundancy  must  be  built  into  the 
system  as  this  type  of  sensor  has  only  marginal  durability  for  ocean- 
ographic  measurements. 
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ABSTRACT 


Thermistor  measurements  of  small  scale  temporal  and  spatial  fluctuations 
of  temperature  were  made  near  the  surface  of  the  sea  from  the  NUC  Oceano- 
graphic  Research  Tower.   Wave  height,  water  particle  velocity,  sound  velocity 
and  salinity  were  measured  simultaneously.   Data  were  subjected  to  statistical 
and  spectral  analyses.   Temperature  fluctuations  were  found  to  vary  with  wave 
height   A  linear  relation  was  found  to  exist  between  temperature  fluctuations 
and  wave  swell.   The  assumed  relationship  was  verified  by  the  results  of 
spectral  analysis.   High  coherence  (.68)  between  wave  height  and  temperature 
fluctuations  occurred  near  the  frequency  of  maximum  wave  energy.   A  relatively 
constant  phase  difference  of  180°  was  observed  in  the  frequency  band  normally 
associated  with  surface  swell  (.03  Hz  to  0.2  Hz).   Water  particle  velocities 
displayed  the  same  relatively  high  coherence  in  the  frequency  band  of  high 
wave  energy.   Temperature  fluctuations  were  180°  out  of  phase  with  the 
u-component  of  velocity  and  90°  out  of  phase  with  the  w-component  of  velocity. 
These  phase  differences  agreed  with  linear  wave  theory.   Coherence  was  observed 
to  be  lower  and  phase  difference  more  random  with  increased  thermal  gradients. 
This  was  attributed  to  an  increase  in  turbulent  temperature  fluctuations. 
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